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Abstract 

The E814-collaboration has found a component of very low m t K + mesons 
with a slope parameter of T ~ 15MeV. We will present a scenario which 
explains the observed slope parameter and which allows us to predict the ex- 
pected slope parameter for kaons produced in heavier systems such as Au+Au. 
Implications for the restoration of chiral symmetry in relativistic heavy ion 
collisions are discussed. 



1 Introduction 



At quark matter '93 J. Stachel presented the first data on kaon spectra obtained by 
the 15-814 spectrometer at the AGS. These data cover only the very lowest range in m t 
(m t — rrik < 25MeV) at somewhat forward rapidities (2.2 < Y < 2.6) and they can be 
parameterized by an exponential in m t with a slope parameter as low as T ~ 15 MeV. This 
slope parameter should be compared with the finding of the E-802 collaboration 0] which 
gives a slope-parameter of ~ 150 MeV at central rapidities for values oim t —m > 100 MeV. 
Also p-Be data, which are essentially equivalent to p-p, show a slope parameter of the 
order of 150 MeV ||. Hence, the slope parameter found by E-814 certainly has to have 
its origin ion the many body dynamics of the heavy-ion collision. 

A steep rise in the spectrum at low transverse momenta has also been observed for 
pions @, |, | . In this case the so called low pt-enhancement can be understood by the 
decay of the delta resonance and event generators which include the deltas such as RQMD 
and ARC || [?J do reproduce the observed spectra very well. In case of the kaons (K + ), 
however, a component due to resonance decay is essentially ruled out, since there very 
few Phi mesons and strange anti-baryons are expected to be produced. 

Another mechanism, which may lead to a cold component in the spectrum is the 



effect of attractive potentials as has been discussed in several papers || || [10], |TT| ■ ^ n 
references M and [TO]] it was shown that the presence of an attractive mean field can affect 
the slope parameter of the particle spectra, in particular the high p t part. If, however, 
the mean field is to affect the soft part of the spectrum, as suggested by ref. |§ for 
the soft pion spectrum at CERN - energy heavy ion collisions, a subtle interplay of the 
expansion time of the fireball and the particles under consideration is needed, as discussed 
in detail in ref. [T^J. This relation essentially is the requirement for adiabaticity of the 
expansion in thermodynamics, meaning the expansion velocity has to be small compared 
to the velocity of the particles to be cooled. While at CERN energies the requirement 
for adiabaticity is not fulfilled, as demonstrated in ref. [Jl~2 , at lower energies such as the 
Bevalac, these conditions hold and a cooling due to the mean field is possible []TD 



In this note we will discuss under which circumstances an adiabatic cooling of the K + 
is possible. We will argue, that the observed very cold component in the kaon spectrum 
requires a rather small expansion velocity of the fireball, which, at the observed tempera- 
tures, can only be achieved if one assumes that the system initially goes through a phase 
with high energy density but low pressure. Such a phase is to be expected from the chiral 
restoration transition. As discussed in a recent analysis of lattice gauge data [13] the chi- 
ral restoration transition involves the decondensation of the gluon condensate which gives 
rise to an effective bag pressure of pbag — — 200MeV/fm 3 . This bag- or decondensation 
pressure lowers the pressure of the system to very small values, nearly zero, and certainly 
much smaller than the 1/3 of the energy density for an ideal gas. As a consequence the 
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system expands much more slowly than an ideal gas, which, for temperature of the order 
of 150 MeV, has an expansion velocity of the order of 1/3 of the velocity of light. At this 
point we should also mention that the lowering of the pressure due to decondensation is 
inherently absent in cascade-type calculations. Those essentially simulate the behavior of 
an ideal gas and thus, as we will show below, predict too large an expansion velocity in 
order for adiabatic cooling to be effective. 

Another ingredient of our model is an attractive potential for the K + . Since the 
K + N scattering amplitude is repulsive simple impulse approximation predicts a slightly 



repulsive optical potential in nuclear matter at zero temperature [14], [15] 



U = -^(1 + ^) a K+N p ~ 25MeV-?- (1) 
2m K M N p 

where a^+N — ~ 0.2 fm is the isospin averaged s-wave kaon nucleon scattering length 
P~6f| . In the framework of chiral perturbation theory []17| , this potential arises from the 



interplay of an attractive scalar piece, which is due to explicit chiral symmetry breaking, 
and a repulsive vector potential as a result of vector meson exchange. If for some reason 
at high temperature the vector meson coupling is reduced, the resulting kaon potential 
could very well turn out to be attractive. Actually lattice gauge results on the behavior 
of the baryon number susceptibility can be understood in a hadronic framework only 
if the vector mesons decouple above T c [f[8| . (For a more detailed discussion see [|19] .) 



It is, however, not the purpose of this article to develop yet another theory about the 
behavior of the K + in matter. We rather want to follow a phenomenological approach 
and determine the strength of the K + optical potential from the available data. As we 
will demonstrate in the following section, this can be done rather unambiguously within 
the framework of our model of collective cooling. Our result for the kaon potential should, 
therefore, be viewed as a phenomenological one, which then needs to be understood from 
underlying theoretical principles. 

This article is organized as follows. In the following section we will present a schematic 
model of the expansion, which explains the essential physics of cooling the kaons and 
shows that the model parameters can be fixed from existing data. In section 3 we present 
results of a more microscopic calculation using relativistic transport theory |2(J and we 



also present out prediction for the larger system Au + Au after having fixed the model 
parameters by comparing with existing data on Si + Au. 

2 A schematic model 

It is widely accepted that in AGS collisions a hot zone consisting mainly of nucleons, 
deltas, and pions is being formed. In our schematic model we assume that these particles, 
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in particular the nucleons and deltas, give rise to an attractive potential for the kaons 
which follows the local baryon density 



U(p(r)) ~ p(r) (2) 

where p(r) is the local baryon density. As a result of the expansion of the fireball, the 
baryon-density and hence the potential change with time. We assume a Woods-Saxon 
type density profile. 

P(r) ~ exp((r-r (t))/A(t)) + l (3) 
so that the potential has the following form 

U(r) = — (4) 

[) exp((r - r (f))/A(f)) + 1 [> 

The expansion of the fireball is then modeled by increasing the radius r and the surface 
thickness A function of time. 

r {t) = r (* = 0) + - 

r 

A(t) = A(t = 0)- ro(t) 



r (t = 0) 

U (t) = U o (t = 0)(^^J (5) 

where v = 1/r is the expansion velocity. By scaling the surface thickness and the cen- 
tral value of the potential with the radius we ensure that the integral over the density 
distribution, i.e. the baryon number remains constant with time. For the Si+Au case we 
assume the following initial conditions for the fireball 



R = 3.5 fm 

A = lfm (6) 

The potential and its time dependence is then determined from eq. (§) using the above 
density function. 

While the motion of the baryons and the resulting mean field potential for the kaons 
is modeled by the time dependence of the density distribution, the propagation of the 
kaons in the potential is treated explicitly by solving the appropriate semiclassical Vlasov 
equation. As it is common practice in the BUU/VUU type (see. e.g. |20|| , transport 
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theoretical description, the Vlasov equation is solved by the so called test particle method, 
in which the phase space distribution is represented by those test particles. The solution 
of the Vlasov equation is then equivalent to solving classical Newton equations of motion 
for these test particles in the potential given above. The initial distribution of these kaon 
test particles follows the same density profile @ in coordinate space. In momentum space 
the kaons are distributed according to a thermal Bose-distribution with temperature 

T=170MeV (7) 

and zero chemical potential. In order to obtain reasonable statistics for the spectra at 
least 50000 test particles are needed. Since we are only interested in the shape of the 
spectra and not the absolute yield, a normalization of the phase space distribution is not 
required. 

Given the above initial condition, we have two parameters left to play with, the ex- 
pansion time r and the depth of the potential Uq. As we will show, both these parameters 
can be fixed by comparison with the data. To this end, let us discuss their effect on the 
spectra. In fig. [I] we show the initial (full histogram) and the final (after expansion) spec- 
trum at fixed potential depth of Uq = —50 MeV and for three different expansion times, 
r = 10 fm/c, r = 5/m/c, and r = 2.5/m/c. We find, that the slope at low transverse 
mass strongly depends on the expansion time used: the longer the expansion time the 
smaller the slope parameter, exactly what one expects from adiabatic expansion. The 
depth of the potential, on the other hand, hardly affects the slopes at low transverse mass 
but rather controls the number of particles contributing to the soft component of the 
spectrum or, equivalently, the value of the transverse mass where the spectrum changes 
to the initial slope. The reason why the depth of the potential hardly affects the slope 
at low transverse masses can be understood classically. Particles are being slowed down 
only by the gradient of the potential, the force, which is effective only at the surface. If 
the system expands too quickly, the surface moves away from the kaons, and thus cannot 
slow them down. If the expansion is slow, on the other hand, only those kaons which 
are trapped inside the potential, will be cooled adiabatically. The number of particles 
'trapped' inside the potential naturally depends on the depth of the potential as only 
those with a kinetic energy smaller than the potential will be kept inside. Therefore, the 
number of particles contributing to the soft component of the spectrum is directly related 
to the depth of the potential. This is demonstrated in fig. [| where we show the resulting 
spectra for the different potentials U Q = -50 MeV, U = -100 MeV, and U = -150 MeV. 
We see that the choice of potential hardly affect the slope of the soft component. 

From the above consideration it is clear, that both model parameters can be fixed from 
an experimental spectrum with a very low temperature component at small transverse 
momenta and a high temperature component at large transverse momenta. The expansion 
time is controlled by the slope of the low temperature component and the depth of the 
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potential by the position at which the spectrum changes from cold to hot. Our schematic 
model, of course, is not realistic enough, to compare with data. To this end, we will use 
a relativistic transport description PDA which takes into account the rescattering of the 



kaons in the fireball as well as the proper transport properties of the particles (nucleon, 
deltas and pions) which make up the fireball. In addition, the transport description makes 
sure, that the expansion of the fireball conserves energy and momentum, which has not 
been addressed in our schematic consideration. 



3 More realistic model 

In this section we want to discuss the transport theoretical calculation of the expansion 
and present the results for the kaon spectra. As already indicated in the introduction, ideal 
gas type expansion is to fast to be effective in cooling the kaons. The important slowing 
down of the expansion is rather a direct result of a transition through a low pressure 
phase which we believe is the chiral restoration transition. This leads to a softening of the 
equation of state and as a result the pressure, which drives the expansion, is lower than 
that of the ideal gas. In order to model the effect of chiral restoration or the softening of 
the equation of state, particles in the fireball should interact via a mean field in addition to 
the standard cascade type collisions. By way of that additional mean field, the equation 
of state can be softened. Thus, in our approach the mean field simulates the effect of 
chiral restoration. 

The model we will use for our calculation is a relativistic transport model (see e.g. 
|]2(J|), which includes nucleon, deltas, pions and kaons. The collisions among these particles 
are controlled by known cross sections, where for the nucleon-delta collisions an improved 



detailed balance method [221 has been used. For the kaon nucleon and kaon delta cross 



section we have assumed a value of 10 mb independent of energy. The K + ir cross section 



we assume to be dominated by the K* intermediate state and we use p3 

00 



&K + 1T 



1 + 4(Vs - m K *) 2 /T 2 K * 



with (T = 60 mb, m K * = 895 MeV, and T = 50 MeV. 

The mean field is based on the Walecka- (a — uj) model, with the possibility to include 
self interaction terms for the scalar field, if needed. This model has the advantage, that 
the necessary field (or decondensation) energy /pressure associated with a change of the 
gluon condensate at the chiral restoration transition |L3) can be simulated by the field 
energy of the scalar sigma-meson at mean field level. It has also been shown, that the 
Walecka model shows a very rapid rise in the energy density to the value of free, massless 



particles at a certain critical temperature, which depends on the parameters used [24 



Of course the Walecka model is not chirally invariant and one may question whether 
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it is the correct tool to study the chiral phase transition. However, since at this point 
we are concerned only about the expansion dynamics, a mean field model which can be 
parameterized such that the pressure of the system is below that of an ideal gas serves 



this purpose. That this is possible in the Walecka model has been shown in ref. |24 . 
Incidentally, using the standard Walecka parameters, the scalar field gives rise to a bag- 
type pressure of ~ 200 MeV/fm 3 , which is about the value extracted in the analysis of 
lattice calculations Hl3| . A more refined calculation, which addresses questions such as a 
temperature dependent vector coupling et c, eventually will have to be done in a chirally 
symmetric model. 

The potential for the kaons is taken of the form 

U K (r) = UA (9) 
Po 

where the density p s is the scalar density < Tpip >. 

Thus, given the energy density as a function of temperature and chemical potential the 
Walecka parameters can be adjusted to reproduce that behavior of the energy density. For 
zero chemical potential fairly reliable and model independent information for the energy 



density can be extracted from lattice gauge calculations [[13|, ^5|, p6fl . However, for HI 
collisions at presently available energies, the energy density needs to be known at finite 
chemical potential. Therefore, lattice results, which are restricted to vanishing chemical 
potential can only be used for guidance and one has to rely on model prediction of the 
properties of the chiral transition at finite chemical potential. Here, however, we want to 
follow the opposite, phenomenological approach. Rather than starting from some model 
input, we will try to extract the properties of the transition from the available data on kaon 
production. As demonstrated in the previous section, it is possible to make statements 
about the expansion velocity and hence of the chiral transition independent of the choice 
of the potential for the kaons, which also can be fixed from the data. With both model 
parameters fixed by data form Si+Au collisions, we then can make a prediction for the 
behavior of the kaon spectra for Au+Au, which have already been taken and are in the 
process to be analyzed p7 |. 



In fig. we show our result for the kaon production data for Si+Au. Together with 
the data, we also show the prediction of the RQMD model determined by the E-814 
collaboration. In order to obtain reliable statistics, in our calculation the kaon spectrum 
has been obtained after averaging over rapidity, while the data are taken at somewhat 
forward rapidities. A detailed analysis of the rapidity dependence will be presented in a 
more extensive article^. Depending on whether or not the last data point is to be taken 

1 In our schematic model (section 2) we have checked that in the limit of a Bjorken type fireball, a soft 
component in the kaon spectrum can be obtained as well, provided the longitudinal expansion is not too 
rapid. Since at AGS energies the fireball formed is somewhere between Bjorken-type and spherical, we a 
rather confident, that in a more extensive simulation, the overall picture will not be changed. 
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into account they are reproduced by a potential of Uo — 50 MeV. The parameters we 
needed for the Walecka model are 

g v = 5.5, g s = 9.27 (10) 

which is to be compared with the original Walecka parameters which have been determined 
from fitting nuclear matter properties 

g° ri9 = 11.56 g°; ia = 9.27 (11) 

Notice that a smaller vector coupling than the original Walecka value is needed in order 
to reproduce the kaon dataQ. This is in qualitative agreement with our arguments for the 
attractive kaon potential, where also a reduced vector coupling is needed. Notice also, that 
a direct comparison with lattice gauge results for the energy density does not determine 
the vector coupling since it does not contribute at zero chemical potential. However, a 
comparison with measured baryon number susceptibilities from the lattice also requires a 
reduction of the vector coupling around T c |L8[], and a more complete calculation will have 



to take the temperature dependence of the vector coupling as given by the lattice gauge 
data into account. In fig. |3| we also show the resulting kaon spectrum obtained without 
any mean field for the baryons. Clearly no soft component has developed in this case. 

Finally, our prediction for the kaon spectrum in case of Au+Au. We simulate the 
Au+Au collision by assuming a fireball of radius 6.5 fm and with 350 baryons inside the 
fireball. In fig. ^ we show the resulting spectrum for this configuration. Again the 
spectrum develops a concave shape but at low m± the slope is not as steep as that for the 
smaller system, mainly because of the increased rescattering of the kaons in the bigger 
fireball. 

Before we conclude we should mention, that in our calculation the protons also show 
a soft component at small transverse momentum with a slope parameter of T ~ 80 MeV 
which is not seen in experiment. We believe that this is an artifact of the simplification 
of our calculation, namely the fact that we have worked with a constant vector coupling. 
In a more realistic calculation, the vector coupling should increase once the system has 
gone through the chiral transition giving rise to additional repulsion in the late stage 
of the reaction. ( The g v = 5.5 we have used is actually an average between the high 
temperature g v — and the low temperature Walecka value ^g^ alecka ~ 11) ) This should 
give the nucleons an additional collective push which will eliminate the soft component 
in the proton spectra. 



2 There is actually some freedom in the choice of the couplings. Both g s and g v can be changed by 
about 25% without affecting the results, provided that their difference is kept at about the same value. 
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4 Conclusions and Outlook 



In this article we have explained the observed cold component in the K + spectrum at AGS 
energy Si + Au collisions by the E-814 collaboration by collective cooling. For the cooling 
mechanism to be successful we had to make three model assumptions. First, that the kaon 
feels an attractive potential in nuclear matter at temperature close to T c . Second, that 
the expansion of the fireball created in the collision is slow, because the system has to go 
through the chiral phase transition. And third, that above T c the coupling of the vector 
mesons is reduced. This reduced vector coupling would then also qualitatively explain 
the attractive kaon potential. 

We have simulated the effect of the chiral transition on the expansion dynamics with 
a relativistic mean field model. In order to slow down the system sufficiently, we had 
to choose a vector coupling which is smaller than that extracted from zero temperature 
properties of nuclear matter. Having fixed the mean field parameters for the nucleons and 
kaons by comparing with the data for Si + Au, we have predicted the K + spectrum for 
the heavier system of Au + Au. Again the spectrum exhibits a concave shape but with a 
larger slope at small transverse mass than that of the small system (Si + Au). Also the 
fraction of soft kaons is somewhat reduced. This is simply a result of the more efficient 
rescattering in the bigger system, which knocks more kaons out of the soft component. 

If the scenario presented in this article is correct, the smaller system of Si+Au seems 
to be better suited to study the effects of chiral symmetry restoration^]. It would, there- 
fore, be desirable if these measurements could be extended over a wider range in p t and 
rapidity, such that the connection with the data of the E-802 collaboration can be made. 
In particular it would be extremely interesting to study the kaon spectra as a function 
of bombarding energy. Within our scenario we would predict that the soft component 
should disappear at lower bombarding energies, because the chiral symmetry restoration 
temperature is not reached and the kaons should not feel an attractive potential. But also 
at higher energies, the soft component eventually should disappear, because for temper- 
atures much larger than T c , when most of the glue is decondensed, the pressure increases 
to that of an ideal gas. As a consequence, the system expands rapidly and the necessary 
condition for adiabaticity is not fulfilled anymore. 

Our present calculation will have to be improved in order to include a dynamical 
change of the vector coupling, which we have taken to be constant. In order to ensure 
energy and momentum conservation, the mean field model has to be modified. This 
modification will be guided by the lattice results on the baryon number susceptibility. 
This also will restrict the choice of mean field parameters and should cure the problem 
with the proton spectra. Furthermore, we are about to investigate, if the measurement of 
HBT correlations of the soft kaons could provide independent experimental information 

3 It may also be possible to create a smaller fireball by triggering on more peripheral Au+Au events. 
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about the long lifetime of the fireball. 

In conclusion, we believe that the soft component in the kaon spectra may very well 
be the first observed signature for chiral restoration in the laboratory. The scenario 
presented in this article may also be applicable to CERN energy reactions, provided that 
the temperature reached is not to far above T c . If this is the case, the analysis of ref. JL2 



has to be repeated and it may very well turn out that the observed soft pions at CERN 

are yet another result of the chiral restoration transition. 
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Figure captions: 

Figure 1: Initial and final kaon spectra for different expansion times and fixed kaon 
potential at U = —50 MeV. For orientation the data and RQMD results of the E-814 
group are also shown. The calculated results are all in arbitrary units and are normalized 
such that the spectra have the same value at m± — m = 0. 



Figure 2: Initial and final kaon spectra for different kaon potentials and fixed expansion 
time r = lOfm/c. For orientation the data and RQMD results of the E-814 group are 
also shown. The calculated results are all in arbitrary units and are normalized such that 
the spectra have the same value at m± — m = 0. 



Figure 3: Kaon spectra from transport calculation. Full line is result including mean 
fields for Baryons and kaons. The dotted line is the result for kaon mean field only, no 
mean field for the baryons and the dashed line is the pure cascade result. The calculated 
results are all in arbitrary units. 
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Figure 4: Kaon spectra for Si+Au (full line) and our prediction for Au+Au (dotted line) 
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